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O ' ABSTRACT 
O 



Context. A new four-telescope interferometric instrument called PIONIER has recently been installed at VLTI. It provides improved 
imaging capabilities together with high precision. 

Aims. We search for low-mass companions around a few bright stars using different strategies, and determine the dynamic range 
p\ ' . currently reachable with PIONIER. 

Methods. Our method is based on the closure phase, which is the most robust interferometric quantity when searching for faint 
companions. We computed fhe,^ 2 goodness of fit for a series of binary star models at different positions and with various flux ratios. 
The resulting x 1 cube was used to identify the best-fit binary model and evaluate its significance, or to determine upper limits on the 
O |i companion flux in case of non detections. 

Results. No companion is found around Fomalhaut, tau Cet and Regulus. The median upper limits at 3cr on the companion flux ratio 
are respectively of 2.3 x 10~ 3 (in 4 hours), 3.5 x 10~ 3 (in 3 hours) and 5.4 x 10~ 3 (in 1.5 hour) on the search region extending from 
5 to 100 mas. Our observations confirm that the previously detected near-infrared excess emissions around Fomalhaut and tau Cet 
are not related to a low-mass companion, and instead come from an extended source such as an exozodiacal disk. In the case of del 
Aqr, in 30 min of observation, we obtain the first direct detection of a previously known companion, at an angular distance of about 
40 mas and with a flux ratio of 2.05 x 10~ 2 ± 0.16 x 10~ 2 . Due to the limited u, v plane coverage, its position can, however, not be 
unambiguously determined. 

Conclusions. After only a few months of operation, PIONIER has already achieved one of the best dynamic ranges world-wide for 
multi-aperture interferometers. A dynamic range up to about 1:500 is demonstrated, but significant improvements are still required to 
reach the ultimate goal of directly detecting hot giant extrasolar planets. 

Key words, techniques: interferometric - binaries: close - stars: low mass - brown dwarfs - planetary systems 

1 . Introduction Adaptive optics on 10-m class telescopes, possibly aided by 

the use of a coronagraph and/or differential imaging techniques, 
Searching for faint companions around bright nearby stars is currently the most wi dely used technique for directly imaging 
is one of the most chal l engin g goals of today ' s astr onomy extrasolar planets (e.g.. Ichauvin et alJ 12004 iKalas et al.l [2008; 
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dOppenheimer & Hinklevl 12009b lAbsil & Mawel l2010h. This iMarois et al.l l2008: Lagrang e*etalJl2009bh . This technique pro- 
/\ ; quest for high-dynamic range observations has been spurred for v ides very high-contrast images at relatively large angular dis- 
^ ! more than a decade by the discovery of more than 500 extrasolar tances from the central star, up to about 1 : 10 6 beyond 1", and can 
planets in our galactic neighbourhood. A proper characterisation detect companions with fl ux ratios around 1 : 100 down to about 
of the physical and atmospheric parameters of the detected plan- 100 mas separations (e.g. iBoccaletti et~aTll2009l) . The upcoming 
ets requires the photon emitted (or reflected) by the planets to generation of high-contrast imagers using extreme adaptive op- 
be disentangled from those emitted by their host star. While a tics on 10-m class telescope promises to further improve these 
temporal separation of the photons is possible in the particular performances. Combined with state-of-the-art coronagraphic de- 
case of transiting systems, most systems require high-dynamic vices such as vortex coronagraphs, the detection of high-contrast 
range imaging with a high angular resolution in order to be char- companions at the diffraction limit (A/D, i.e., abou t 50 mas at K 
acterised. band on a 10-m telescope) will be enabled (e.g.. IMawet et al.l 

l201ll) . 



* Based on observations obtained at the European Southern Pushing the discovery space within the diffraction limit of a 

Observatory (ESO) Very Large Telescope Interferometer (VLTI), single aperture requires interferometric techniques. On a single 

Paranal, Chile. pupil, aperture-masking techniques have recently allowed dy- 

** Postdoctoral Researcher F.R.S.-FNRS (Belgium), email: namic ranges up to 1:1000 to be reached down to the diffrac- 

absil@ast.ro . ulg . ac . be. tion limit, and permitted decent contrasts (~ 1 : 100) to be reached 
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down to A/3D (iHinklev et alJ201 ltlLacour et al.l20Tlh . Another 
example is the Palomar Fiber Nuller, which also works on a sin- 
gle telescope, and has reached a 1:500 dynamic range around 
the bright star Vega in a search region extending significantly 
within the diffraction li mit of the Palomar Hale 5-m telescope 
dMennesson et alj201 ll) . To extend the search region even closer 
to the central star, interferometry with multiple apertures sepa- 
rated by more than 10 m is mandatory. Interferometry on mul- 
tiple apertures has, however, not yet reached the same dynamic 
range as adaptive optics or even aperture-masking ob servations, 
at lea st when the target star is mostly unresolved dAbsil et al.l 
120101) . The dynamic range can be significantly increased by 
obser ving resolved targets (e.g., iDuvert et alJl2010l IZhao et all 
1201 ll) . but such targets are scarce in the context of main sequence 
stars observed on 100-m class baselines. 

Here, we present the first results of a new 4-telescope in- 
terferometric instrument called PIONIER (Precision Integrated- 
Optics Near-infrared Imaging ExpeRiment), which has recently 
been installed and commissioned at th e Very Large Telescope 
Interferometer (Le Bou quin et al.ll20lTb - PIONIER aims at pro- 
viding enhanced imaging capabilities at the VLTI, together with 
an improved accuracy and stability of interferometric observ- 
ables (visibility and closure phase). One of its main scientific 
goals is to directly detect faint companions around bright stars 
(up to a magnitude H ~ 7). The purpose of this paper is to dis- 
cuss the dynamic range of this new instrument based on the ob- 
servations obtained during the commissioning phase and the first 
scientific observing runs. 

2. Observations and data reduction 

The observations presented in this paper were obtained during 
the commissioning phase and the first scientific observing runs 
of PIONIER at VLTI, which took place in Fall 2010. A large 
number of stars has been observed during several nights from 
October 25 to December 25, using only the 1.8-m Auxiliary 
Telescopes of the VLTI. Within this sample, we selected four 
stars that are representative of two possible observing strategy: 
either a long integration for a deep search for companions (cases 
of Fomalhaut and tau Cet) or a short ("snapshot") integration for 
a shallow survey-type search (cases of del Aqr and Regulus). The 
observing log related to these four stars is displayed in Table Q] 
Observations of the science targets were interleaved with obser- 
vations of reference stars to calibrate the instrumental contribu- 
tion in the observed quantities. Calibrators were usually chosen 
close to the science targets, both in terms of position and magni- 
tude. The main parameters (magnitudes, diameters) of the four 
targets and their calibrators are listed in Table [2] Because some 
observations were obtained during commissioning, nearby cal- 
ibrators were not always available. In these cases (e.g., on the 
night of October 28), we carefully checked on long calibration 
sequences that the instrumental transfer function does not de- 
pend on the target position, which was always the case. 

Various baselines were used during the observations, with 
ground separations ranging from 16 m (D0-E0, E0-G0) to 128 m 
(A0-K0). The mean ground baseline lengths of the four con- 
figurations used during the reported observations are of 44 m 
(E0-G0-H0-I1), 53 m (D0-E0-H0-I1), 63 m (D0-G1-H0-I1), and 
86 m (A0-G1-K0-I1), respectively. The observations are divided 
into observing blocks (OBs), corresponding to the lines of 
Table [1] Each OB is composed of several files, which con- 
sist in 100 successive fringe measurements and typically last 
1 min. The total amount of time spent on each of the four 
targets, including calibration observations, amounts to about 



Table 1. Observing log for Fomalhaut, tau Cet, del Aqr, and 
Regulus. 



Night 


MJD 


Files 


Star 


Baseline 


2010 Nov 30 


55530.999 


3 


HD 209688 


E0-G0-H0-I1 




55531.007 


10 


Fomalhaut 






55531.023 


5 


HD 209688 






55531.035 


10 


Fomalhaut 






55531.048 


5 


HD 219784 






55531.055 


10 


Fomalhaut 






55531.067 


10 


HD 219784 






55531.098 


10 


Fomalhaut 






55531.114 


10 


HD 219784 




2010 Dec 6 


55537.047 


5 


HD 219784 


D0-G1-H0-I1 




55537.058 


5 


Fomalhaut 






55537.067 


5 


HD 219784 






55537.076 


5 


Fomalhaut 






55537.084 


5 


HD 219784 






55537.090 


5 


Fomalhaut 






55537.096 


5 


HD 219784 




2010 Dec 3 


55534.053 


5 


tau Cet 


E0-G0-H0-I1 




55534.060 


5 


HD 8705 






55534.066 


5 


tau Cet 






55534.075 


5 


HD 14376 






55534.081 


5 


tau Cet 






55534.092 


5 


HD 15694 






55534.106 


5 


HD 16526 






55534.113 


5 


tau Cet 






55534.121 


5 


HD 15694 




2010 Dec 19 


55550.031 


5 


HD 8705 


A0-G1-K0-I1 


... 


55550.041 


5 


tau Cet 




... 


55550.051 


5 


HD 14376 






55550.073 


5 


tau Cet 






55550.083 


5 


HD 8705 




2010 Oct 28 


55498.089 


3 


HD 221745 


D0-E0-H0-I1 




55498.096 


3 


del Aqr 






55498.105 


3 


gam Tuc 




2010 Nov 28 


55529.319 


4 


HD 45018 


E0-G0-H0-I1 




55529.335 


5 


Regulus 






55529.343 


3 


HD 83425 






55529.351 


3 


Regulus 






55529.359 


3 


HD 82381 






55529.369 


3 


Regulus 





4h (Fomalhaut), 3h (tau Cet), 30 min (del Aqr), and lh30 
(Regulus), respectively, for a total number of calibrated OBs of 
7 (Fomalhaut), 6 (tau Cet), 1 (del Aqr) and 3 (Regulus). The ob- 
serving conditions generally ranged from good (seeing of 0'.'8, 
to = 5 msec) to fair (seeing of 1'.'2, To = 2 msec), with the ex- 
ception of the observations of 2010/1 1/28 on Regulus, where the 
conditions were poor (seeing of 1'.'4, to = 1 .5 msec). 

The interferometric concept and the data reduction strat- 
egy of PIONI ER are directl y inspired from the IONIC-3 ex- 
periment (Ber ger et al.l 12001). The interested reader can find a 
detailed description of the instrument and its performance in 
iLe Bouquin et al.l d201 lb . In brief, PIONIER relies on a pairwise 
co-axial combination embedded inside an integrated optics com- 
ponent. Before the flux is injected into the chip, the four beams 
are modulated at non redundant velocities (— 3v, -v, +v, +3v) to 
generate the six fringe signals. The length of the shortest scan 
is generally set to 80 yum. For bright stars, v is imposed by the 
detector frame rate and the need to correctly sample the fringe 
at the highest frequency f„ mx = 6v/A. In the observations pre- 
sented here, the velocity was set to v = 70^m/s, which corre- 
sponds to fringes at 80, 160, and 240 Hz. This is just fast enough 
to freeze the turbulent optical path delay (OPD) under median 
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Table 2. Main parameters of our four targets and their calibra- 
tors. 



Star 


Type 


H 


Old 


Ref. 


Fomalhaut 


A3V 


1.03 


2.223 ± 0.022 


1 


HD 209688 


K3III 


1.58 


2.71 ±0.030 


2 


HD 219784 


G8III 


2.01 


2.13 ± 0.025 


2 


tau Cet 


G8V 


1.72 


2.015 ± 0.011 


3 


HD 8705 


K2III 


2.17 


2.06 ± 0.040 


2 


HD 14376 


K5III 


2.92 


1.457 ± 0.019 


4 


HD 15t>94 


V 1 TTT 


z.ol 


1. /DO ± U.Uzj 


4 


HD 16526 


K4III 


3.12 


1.389 ±0.017 


4 


del Aqr 


A3V 


3.14 


0.825 ± 0.012 


5" 


HD 221745 


K4III 


3.06 


1.426 ±0.019 


4 


gam Tuc 


F4V 


3.00 


0.998 ± 0.070 


6 


Regulus 


B7V 


1.57 


1.47 ± 0.20 


7 b 


HD 45018 


K5III 


2.00 


2.50 ± 0.041 


2 


HD 83425 


K3III 


1.68 


2.64 ± 0.028 


2 


HD 82381 


K2III 


2.07 


2.16 ±0.026 


2 



No tes. Reference s : (1) [Absi et all j2009t> . (2) iBorde et all f2 002). 
(31 iDi Folco etai] d2007h i . (4) iMerand etaU d2005h. ( 5) iKervella etal] 
d2004h . (6) lLafrasse et alJd2010h . (7) lMcAlister et aT]d2005h . 
(a) Diameter based on surface-brightness rel ationships, using m agni- 
tudes V = 3.269 ± 0.01 and K = 3.078 ± 0.02 dBouchet et aD 19911) . 
(t) The uncertainty takes into account the stellar oblateness. 

atmospheric conditions. The outputs of the chip are dispersed 
over seven spectral channels across the H-band, providing low 
{R = 35) spectral resolution. During the observation, each scan is 
processed by a quick-look algorithm that also implements slow 
group tracking. 

All th e data have been reduced with the pndrs package de- 
scribed in lLe Bouquin et all (1201 ll) . For each file and each spec- 
tral channel, the pndrs package provides six visibility and four 
closure phase measurements. The statistical uncertainties are es- 
timated by the dispersion over the 100 scans contained in each 
file, and typically range from 0.25 to 3 degrees for the closure 
phase, depending on the target brightness and atmospheric con- 
ditions. The statistical uncertainty includes contributions from 
the classical detection noises (detector and photon noise), as well 
as from the atmospheric OPD fluctuations faster than the scan 
rates. The instrumental and atmospheric contributions to the vis- 
ibilities and closure phases (the so-called transfer function) is 
monitored by interleaving the observation of science stars with 
calibration stars. Some nights show an extremely stable transfer 
function, down to 0.1 deg, while other nights require fitting the 
transfer function by a drift, indicating that we are facing some 
non-stationary biases. This aspect requires additional commis- 
sioning before being fully understood. 

3. Data analysis strategy 

In this section, we discuss how the PIONIER calibrated data 
set is used to search for faint companions, compute confidence 
levels for tentative detections, and derive upper limits to the 
presence of off-axis companions in case of non-detections. Our 
analysis is based solely on the calibrated closure phases, which 
have two significant advantages over (differential) visibilities 
and (differential) phases. On the one hand, at first order, clo- 
sure phases are not affected by phase dis turbances ahead of 
the beam combiner (see e.g., Monnier 2003). Therefore, closure 
phases are not affected by atmospheric turbulence and other in- 
strumental effects such as mechanical vibrations, which strongly 
reduce the accuracy of other observable quantities. On the other 



hand, closure phase are especially sensitive to faint off-axis com- 
panions, as they differ from zero only for non point-symmetric 
targets. One of the associated benefits is that the detection of 
a faint companion does not depend on the knowledge of the 
parent star's photospheric parameters, provided that it is point- 
symmetric; however, poor knowledge of the primary stellar di- 
ameter may lead to larger uncertainty on the photometry of a 
potential companion. 

Our companion search strategy is based on the computa- 
tion of the x 2 distance between our data set and several mod- 
els where a faint companion is added to a user-specified pho- 
tospheric model for the central star. The x 1 goodness-of-fit is 
computed for a series of angular separations along east and north 
between the primary star and its companion. The first step in this 
analysis is therefore to define the region around the star in which 
companions will be searched for. 

3. 1 . Determining the search region 

There are three limitations to the maximum angular separation 
at which a companion can be searched for in the PIONIER data. 
The first one is related to the use of single-mode fibres in the 
instrument. The coupling efficiency of point-like sources into a 
single-mode fibre depends on the angular separation between the 
source and the optical axis. The transmission profile, which re- 
sults from the overlap integral between the turbulent image of 
a point-like source and the fundamental mode of the fibre can 
be approximated by a Gaussian profile, whose full width at half 
maximum (FWHM) depends on the fibre parameters, the tele- 
scope diameter, the observing wavelength, and the strength of 
the atmospheric turbulence. In the case of PIONIER, working 
with the 1.8 m ATs under typical atmospheric conditions (see- 
ing of 078 in the visible), the FWHM of the Gaussian profile is 
about 400 mas. The transmitted flux of off-axis companions lo- 
cated at angular separations larger than 200 mas from the central 
will therefore be attenuated by more than 50%. 

A second limitation to the angular distance of detectable 
companions comes from the separation of the fringe packets as- 
sociated to the two stars in terms of optical path delay (OPD). 
To properly detect a potential companion, we ask the fringe 
packet separation to be less than half the size of the smallest 
scan, which amo unts to 40 /zm for the ba selines with the slowest 
scanning speed dLe Bouquin et al.ll201 lh . The OPD separation 
of the fringe packets is directly related to the angular separa- 
tion AO of the two objects in the sky: AOPD = B AO cos 0, with 
B cos the projected baseline length. Using a mean projected 
baseline of about 40 m for the compact configuration (E0-G0- 
H0-I1, used for most targets) and taking AOPD max = 20/im, 
we end up with A0 mdx - 100 mas. A fringe packet separation 
smaller than 20 fim also ensures that the two fringe packets are 
(at least partially) superposed in all cases, as the size of the fringe 
packet is A 2 /AA 63//m when the H-band signal is dispersed 
onto seven spectral channels. 

The last limitation to the size of the search region comes 
from the spectral sampling of the closure phase signal. Because 
the angular resolution of the array is proportional to the observ- 
ing wavelength, the closure phase signature of an off-axis com- 
panion oscillates around zero as a function of wavelength. These 
oscillations must be properly sampled by the spectral resolution 
of the instrument in order to produce a uniq ue (non-alia s ed) so - 
lution in the fitting process. As discussed bv lAbsil et al.l (12010). 
the period in the closure phase signal, which is roughly given 
by P,\ = A 2 /(B AO - A), must be larger than four times the spec- 
tral channel size. In the present case, considering a mean sky- 
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projected baseline B — 40 m for the compact configuration (the 
most frequently used in our observations) and taking the chan- 
nel width of 0.045 fim in the PIONIER dispersed data into ac- 
count, this constraint translates into A# max =; 72 mas. This does 
not prevent detecting companions located further away, but these 
companions can possibly produce aliases within the inner 72 mas 
region and therefore lead to an ambiguous position. 

Taking these three limitations into account, we restrict our 
search region to a radius of about 100 mas around the target stars. 
Detections within this region should be mostly robust and unam- 
biguous. It must, however, be noted that companions located fur- 
ther away (within about 500 mas from the central star) could lead 
to "false" detections within our 100 mas search region. The pres- 
ence of such companions would, however, generally be known 
from AO-assisted (coronagraphic) observations on 10-m class 
telescopes (see Sect. [TJ. 



values of the first parameter around the minimum. The result- 
ing^ is then translated into a probability, and the error bars are 
defined as the 68% confidence interval around the minimum for 
each individual parameter. 

Finally, we convert the 2D map of the minimum renor- 
malised x 2 into a 2D probability map, following the same pro- 
cedure as in Eq.Q] 



P™ P (x,y)= 1 -CDF V 



ymin,.(y 2 (x,v,r) 



(2) 



A visual inspection of this map is performed to check for mul- 
tiple solutions. We also inspect the closure phase data and the 
best-fit closure phases to look for possible outliers or systemat- 
ics effects on our data set, which could produce artefacts in our 
search results. 



3.2. Searching for companions 

Our search strategy consists in building a series of binary models 
with off-axis companions at all possible locations (x,y) within 
the search region and with various flux ratios r, starting at 
r — to include the single-star model. Our binary models fea- 
ture a realistic description of the primary stellar photosphere, 
based on interferometric measur ements where avail able or on 
surface-brightness relationships dKervella et alj |2004). They in- 
clude the effects of limb-darkening and of a possible elongation 
of the photospher e due to rapid stellar rotation, as described in 
Absil et al] (120081) . The off-axis companions are supposed to be 
fully unresolved by the interferometer (i.e., point-like) and are 
thus only parameterised by their flux ratio r with respect to the 
primary star. For each individual model (x, y, r), we compute 
the associated closure phases at the relevant observing dates and 
wavelengths, and derive the^ 2 goodness-of-fit to the whole data 
set. This quantity is divided by the number of degrees of free- 
dom v of our^ 2 variable to obtain a cube of reduced^ 2 , noted 
X 2 (x,y,r). The search strategy is then based fully on this cube 
and on the assumption that our data set follows purely Gaussian 
statistics. A 2D illustration of the^ 2 cube for our four targets is 
given in Fig.Q] where the minimum^ 2 as a function of flux ratio 
has been selected at each position. In these maps, the maximum 
value in the scale always corresponds to the single-star model. 

The following step is to search for the global minimum of 
the xl cube and to renormalise the cube so that ^ 2 = 1 for the 
global best-fit binary model. In doing so, we implicitly assume 
that this best-fit model actually corresponds to the true nature of 
the observed source. Based on the renormalised;^ 2 , we compute 
the probability Pq for a random variable in the x 2 r distribution 
to be equal to or larger than the x 2 r associated to the single-star 
model (r = 0): 



P = 1 - CDF V 



vxfy = 0) 



r(X r (x,y,r)) 



(1) 



where CDF V is the x 2 cumulative probability distribution func- 
tion with v degrees of freedom. If Pq is below a pre-defined 
threshold, we can reject the single-star model and the best-fit bi- 
nary solution is considered as significant. The threshold is gener- 
ally fixed at a 3<x level, i.e., at a probability of 0.27%. The error 
bars on the best-fit flux ratio and position are then determined 
by inspecting the dependence of the x 2 around the global min- 
imum with respect to each of the three parameters. In practice, 
we fix one parameter, search the minimum^ 2 with respect to the 
two other parameters, and repeat the same procedure for several 



3.3. Deriving sensitivity limits 

In case no companion is detected, we derive sensitivity limits as 
follows. The^ 2 cube is first renormalised so that^ 2 = 1 for the 
single-star model. This corresponds to assuming that there is ac- 
tually no companion around the observed star. The renormalised 
X 2 cube is then converted into a probability cube: 



P m Ux,y,r) = l-CDF y 



vX 2 (x,y,r) 
X 2 r (r = 0)) 



(3) 



For each position (x, y), we start at r — and increase the flux 
ratio of the off-axis companion until the associated probability 
becomes lower than a pre-defined threshold, usually 0.27% for 
a 3<x significance level. If such a companion had been present in 
the data, the single-star model would have been inconsistent with 
the data at 3cr or more, and a detection would have been reported. 
This flux ratio is the lowest that could have been detected at that 
particular position. By repeating this procedure at all positions, 
we create a sensitivity map giving the 3<x upper limit on the flux 
ratio of potential companions: 



S mdp {x,y) - r , so that P cu be(x,y, r) = 0.27% 



(4) 



This map can then be used to produce (cumulated) histograms 
of the sensitivity on the search region, and define the sensitivity 
limit for various completeness levels. 

4. Results 

4.1. Long integrations on Fomalhaut and tau Cet 

Deep integrations, of about 3h and 4h respectively (including 
calibration time), have been obtained on Fomalhaut and tau Cet. 
We used the whole data set for each star to compute the mini- 
mum^ 2 maps, which are displayed in Fig.Q] before any renor- 
malisation. The number of degrees of freedom for the^ 2 variable 
is respectively 1 154 and 629. With such large data sets, it is ex- 
pected that any potential systematic error(s) are mostly averaged 
out, leaving a clean Gaussian statistics. Indeed, it must be noted 
that thex 2 values are close to unity in both cases, which confirms 
that error bars have been properly evaluated by the data reduc- 
tion software and which suggests that systematic errors do not 
dominate the overall error budget. A quick inspection of the best 
X 2 maps shows a relatively high uniformity of the best x 2 (nar- 
row range of values), and the absence of particularly outstanding 
minima. 

To properly assess the presence of significant minima in the 
maps, we renormalise the x 2 cube so that x 2 — 1 for the best-fit 
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100 50 -50 -100 100 50 -50 -100 

RA [mas] RA [mas] 



Fig. 1. Two-dimensional representation of thex 2 cubes, showing the minimum x 2 obtained at each point of the search region when 
exploring the flux ratio of our binary star models in an attempt to reproduce the Fomalhaut (top left), tau Cet (top right), del Aqr 
(bottom left), and Regulus (bottom right) data sets. These maps are computed before any renormalisation of the;^ cube. 



model, as described in Section [3~2l The single-star model is then 
associated with ax 2 of 1.05 for Fomalhaut and 1.14 for tau Cet. 
The corresponding probability for the single-star model to repro- 
duce the data is 10% (or 1 .6cr) and 0.7% (or 2.1 o~), respectively. 
The minima of the x 2 maps can therefore not be considered as 
significant, and we conclude that no firm detection is reported 
around these two stars. 

Following Eq.|U we derive the 3cr upper limits on the flux of 
off-axis companions within the search region around both stars. 
The resulting sensitivity maps are illustrated in Fig. [2] (top and 
middle left). They show sensitivity limits in the range 0.10%- 
0.34% for Fomaihaut and 0.1 1%-0.57% in the case of tau Cet, 
suggesting that companions as faint as 1:1000 of the stellar 
flux would be detectable at some particular positions within the 
search region. The cumulated histogram of the sensitivity level 
across the whole search region is illustrated by the solid curves in 
the left-hand side plots of Fig. [5] The median and 90% percentile 
sensitivity limits are 2.3 x 10~ 3 and 2.8 x 10~ 3 (Fomalhaut) and 
3.5 x 10~ 3 and 4.5 x 10~ 3 (tau Cet), respectively. Another, more 
conventional way to represent the sensitivity of our observations 
is to consider concentric annuli, in which cumulated sensitivity 
histograms can be inspected to search for, e.g., median and 90% 
completeness levels. This is illustrated on the right-hand side of 
Fig. [2] (top and middle). The full dynamic range of PIONIER 



is achieved as close as 5 mas for Fomalhaut and 3 mas for tau 
Cet, and is maintained all the way to the edge of the search zone 
(100 mas), where the effect of the Gaussian profile of the fibre 
transmission starts to be noticeable. 



4.2. Short integrations on del Aqr and Regulus 

To illustrate the capabilities of short integrations ("snapshots") 
in terms of companion detection limits, we have selected two 
representative targets observed during the first scientific runs of 
PIONIER: del Aqr and Regulus. The former was observed on 
only one occasion (see Table H}. The data set shows very obvi- 
ous non-zero closure phases (see Fig. |3). The associated min- 
imum x 2 ma P (Fig- H bottom left) shows that the single-star 
model is inconsistent with the data set, giving^ = 8.2 before 
any renormalisation of the error bars. This x 2 increases to 14 
once the error bars have been renormalised to obtain x 2 = 1 for 
the best-fit binary model. With 53 degrees of freedom in the^ 2 
distribution, the "null hypothesis" (absence of companion) can 
then be rejected at about 40cr. In Fig. [4] we show the proba- 
bility map defined in Eq. [2] Owing to the very incomplete u, v 
plane coverage provided by a single OB (see inset in Fig. |4), 
the position of the detected off-axis object cannot be unambigu- 
ously determined. Three positions are almost as likely (mini- 
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Fig. 2. Left. Map of the 3<x upper limit on the flux ratio of companions around Fomalhaut (top), tau Cet (middle) and Regulus 
(bottom). The uppermost 1% values have been clipped to reduce the colour scale range. Right. Associated sensitivity as a function of 
angular distance, for two completeness levels (50% or 90%). The 50% completeness level corresponds to the median 3cr sensitivity. 



mum^ ranging from 1.00 to 1.16, i.e., at less than l<x from 
each other), and are respectively located at (E: 31.0 + 0.10 mas, 
N: 21.0 ± 0.26 mas), (E: 28.6 ± 0.18mas, N: 29.4 + 0.34 mas), 
and (E: 26.6 + 0.12 mas, N: 38.2 + 0.27 mas). All other posi- 
tions have probabilities below 10~ 3 . The best-fit flux ratio, which 
amounts to 2.05 x 10~ 2 + 0.16 x 10 -2 , is compatible within lcr 



for the three possible positions. The closure phases associated to 
the three best-fit solutions are displayed in Fig. [3] 

In the case of Regulus, three successive OBs have been ob- 
tained (see Table[TJ. This is representative of the quantity of data 
one would get in a large survey for companions around main 
sequence stars, although the number of files per OB (3 to 5) 
was rather low in the present case. The analysis of this data set 
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Fig. 3. Calibrated closure phases for the four triangles and the 
three MJDs of the del Aqr data set. The closure phases of the 
three best-fit binary models are represented by solid, dotted, and 
dashed lines. 



with the x 1 method results in the absence of significant detec- 
tion: once the x 2 cube has been renormalised as discussed in 
Section [3721 the single-star model has xl = 1.14, which cor- 
responds to a probability of 7.6% (equivalent to 1.8<x) to re- 
produce the data set, taking the 230 degrees of freedom into 
account. The single-star model can therefore not be rejected. 
We compute sensitivity limits following Section 13.31 leading to 
the sensitivity map displayed in Fig. [2] (bottom left). The bright 
lines within the map are due to the limited u, v coverage, which 
leads to low sensitivities at some particular positions within the 
search region. The median and 90% percentile sensitivity lim- 
its computed across the whole search region are respectively of 
5.4 x 10" 3 and 7.9 x 10~ 3 . The corresponding sensitivity curves 
with respect to angular separation are shown on the right-hand 
side of Fig.|2](bottom), for 50% (median) and 90% completeness 
levels. 
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Fig. 4. Probability (Eq. [2]i for the best-fit binary models to re- 
produce the data set on a 60 x 60 mas search region around del 
Aqr. The highest three maxima are identified by circles. The 
u, v coverage associated to this data set is represented in the up- 
per right inset, including wavelength dependence (using colours 
from blue to red), with the axes graduated in cycles per arcsec. 

5. Discussion 

5.1. Double-blind search for fake companions 

To validate the sensitivity limits computed in Section l4~Tl for the 
deep integrations, we performed a double-blind test to search 
for fake companions in the data sets where no companion was 
initially detected. One of the co-authors introduced into the cal- 
ibrated data set a series of fake companions of various flux ra- 
tios, located at various positions within the search region. The 
range of angular separations was limited to a minimum of 5 mas 
to ensure an optimum sensitivity, and to a maximum of 50 mas 
to reduce computation time (extending it to 100 mas would not 
have significantly changed the results). Then, the first author, 
unaware of the flux ratios and positions of the fake companions, 
performed a search for companions on the new data sets using 
the procedure described above. In all cases, the overall best-fit 
binary model was identified and the x 2 r cube renormalised so 
that this best-fit model corresponds to x 2 - 1 ■ The probabil- 
ity that the single-star model would reproduce the data set was 
then computed and converted into a number of sigmas, as in 
Section 13.21 The flux ratio and position of all best-fit models 
were finally compared with the input fake companion parame- 
ters. 



5.1 .1 . Long integrations on Fomalhaut and tau Cet 

The blind tests were arranged in series of ten test cases for sev- 
eral values of the flux ratio. A total of 21 different flux ratios 
ranging between 0.05% and 0.6% were investigated in the case 
of Fomalhaut, while seven different flux ratios between 0.15% 
and 0.6% were used in the case of tau Cet. The results of the 
blind tests are illustrated in Fig. [5] The left-hand side plots show 
the fraction of companions that was correctly identified (position 
and flux ratio) as a function of the flux ratio (i.e., the detection 
"completeness"). The agreement with the x 2 analysis is excellent 
in the case of tau Cet, while the blind test suggests a sensitivity 
slightly better than expected from the x 2 analysis in the case of 
Fomalhaut. In the latter case, the median sensitivity within the 
5-50 mas region would be around 0. 19% instead of the expected 
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0.23%. The right-hand side plots display the measured flux ra- 
tios of the best-fit binary models as a function of the true flux 
ratios. These plots demonstrate the absence of bias in the deter- 
mination of the companion brightness. They also confirm that 
the flux ratio of false detections generally does not exceed 0.2% 
and 0.35%, respectively, for Fomalhaut and tau Cet. 

5.1.2. Short integration on Regulus 

A similar analysis was performed for Regulus, with contrasts 
ranging from 0.1% to 2.0% for the fake companions. Because 
the m, v plane coverage is less dense in this case, several replica 
of the fake companions introduced in the data sets appear fre- 
quently in the x 2 cubes, at various positions in the search re- 
gion. These can be seen as the side lobes of the instrumental 
PSF achieved during this snapshot. Due to the statistical noise 
and systematic errors (which are not fully averaged out due to 
the low number of OBs/files), the most significant minimum in 
the x 2 cube does not always correspond to the actual position 
of the companion. In Fig. O we have therefore used two defini- 
tions of the completeness level: a "detection" is reported when 
a significant minimum (i.e., at more that 3<x) is observed in the 
X 2 cube at the companion position, while the "position" of the 
companion is deemed found only when the global minimum in 
the cube is located less than 1 mas away from the actual com- 
panion position. The left-hand side plot of Fig. [6] illustrates a 
behaviour that was already noted in the case of del Aqr: even 
when the detection is (very) clear, the companion position gen- 
erally remains ambiguous due to the limited u, v plane coverage, 
and the global best fit is not always found at the right position. 
This is also the reason for the presence of filled red dots in the 
right-hand side plot for contrasts ranging from 0.5% to 0.9%: the 
estimated companion position is generally wrong in this contrast 
range. Nonetheless, the results of the blind test in terms of pure 
detection (disregarding position) are fully compatible with the 
X 2 analysis. 

Overall, we have thus used three different approaches to eval- 
uate the sensitivity levels, which all give similar results: 

- the x 2 analysis of the original data set, using a 3cr detection 
limit, 

- the detection rate (completeness) in the blind tests, 

- the typical flux ratio for the non-significant detections (< 3<x) 
in the blind tests. 

The agreement between the various approaches further validates 
the sensitivity levels and confirms that 3cr is a reasonable signif- 
icance level for identifying candidate companions. It also vali- 
dates a posteriori our assumption of Gaussian statistics for our 
data sets. 

5.2. Notes on individual targets 

Fomalhaut Fomalhaut is a bright A4V star located at 7.7 pc. 
Near-infrared interferometric observations with VLTI/VINCI 
have recently revealed a K-band excess emission close to the 
photosphere with a flux ratio of 8.8 x 10~ 3 + 1.2 x 10~ 3 
d Absil et al.ll2009l) . In that paper, it was argued that the excess 
emission most probably comes from an extended source, rather 
than from a point-like companion. The absence of closure phase 
signal in our PIONIER data further confirms this conclusion, and 
allows the presence of any companion with a flux ratio greater 
than 2.8 x 10~ 3 to be rejected at 3<x with a 90% completeness 
level on a 100 mas field-of-view. According to the models of 



iBaraffe et al.l (Il998l) . this corre sponds to a mass li mit around 
0.175M o for an age of 290 Myr (IDi Folco et al.l2004 . We there- 
fore confirm that the previously detected excess emission at K 
band comes from an extended (and mostly point-symmetric) 
source. 



tau Cet This G8 main sequence star located at 3.6 pc is also 
known to have a K-band excess emission, with a flux ratio of 
9.8 x 10~ 3 + 2.1 x 10^ 3 derived from high-precision visibility 
measu rements at the CHA RA array with the FLUOR instru- 
ment (Di Folco et al. I l2007h . As in the case of Fomalhaut, our 
PIONIER observations confirm the absence of low-mass com- 
panion within the first 100 mas around tau Cet, with a 90% com- 
pl eteness level of 4.5 x 10 3 at 3<x. According to the models 
of Barafi Vet al.l (Il998l) . this corresponds to a ma ss limit around 
0.09M o for an age of 10 Gyr (IDi Folco et al.ll2004l) : however, the 
significance of the best-fit binary model for the tau Cet data set is 
rather close to 3cr, so that the possibility that a faint companion 
is actually present cannot be completely ruled out. The best-fit 
companion would have a flux ratio of 3.0 x 10~ 3 and would be 
located at (E: 16.9 mas, N: 49.5 mas) from tau Cet. This poten- 
tial companion would, how ever, not be bright eno ugh to explain 
the K-band excess found bv lDi Folco etall ( l2007l) . 

del Aqr This A3 main sequence star, located at 49 + 4pc, is a 
know n astrometric binary with poorly constrained orbital param- 
eters dGoldin & Makarovll2007l) . lLagrange etaf] (l2009al) shows 
that its radial velocity is variable and also identifies it as a bi- 
nary, but fails to constrain the orbit of the companion due to 
the small time span of their observations (< 400 days). With 
our PIONIER observations, we directly detect the companion 
for the first time, although with some ambiguity on its position. 
The orbital parameters cannot be refined based on this sole mea- 
surement, and would require the binary to be observed again at 
several phases along its orbit. However, with this single snap- 
shot, we can readily estimate the spectral type of the compan- 
ion. The H-band flux ratio amounts to 2.05 x 10~ 2 + 0.16 x 10~ 2 , 
which corresponds to A H _ 4.22 + - 09. Ta king the magni- 
tude (H = 3.14 + 0.02, iBouchet et al.|[l99Tl) and distance of 
del Aqr into account, the companion has an absolute magnitude 
Mh — 3.89 + 0.20. Assuming that the detected point-like source 
is a main-sequence star orbiting del Aqr, the companion would 
then have a spectral type around G5V. 

In an attempt to lift the ambiguity on the companion posi- 
tion, we have tried to fit the wavelength-differential squared vis- 
ibilities together with the closure phases, using the exact same 
method. Based on our experience with PIONIER data, the differ- 
ential visibilities are generally less constraining than the closure 
phases when searching for faint companions, for the reasons ex- 
plained in Sect. [3] However, the additional constraints brought 
by the differential visibilities seem to identify the middle posi- 
tion in Fig. [4] as the true solution. The angular separation and 
position angle would then be 4 1 .02 + 0.37 mas and 44.2 + 0.5 de- 
grees, respectively. The ratio of probabilities between the three 
possible solutions is, however, not high enough to give a defini- 
tive answer. 



Regulus Our observations rule out the presence of companions 
with flux ratio greater than 7.9 x 10 -3 at 3cr with a 90% com- 
pleteness level on a 100 mas field-of-view around this B7 main 
sequence star. According to the models of Bara ffe et al.l ([1998), 
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Fig. 5. Blind test result for Fomalhaut (top) and tau Cet (bottom). Left. Cumulated histogram of the 3cr detection limit for companions 
located within the 100 mas search region, based on the^- 2 analysis. Over-plotted in diamonds are the results of the double-blind 
tests discussed in Section IBTTl Right. Measured vs. true flux ratio for the blind tests. Filled points are used when 50% or more of 
the companions have been found at the good position and with the good flux ratio (empty points otherwise). Green points are used 
when the mean significance of the best fit is larger than 3<x (red points otherwise). The error bars represent the dispersion of the 
best-fit flux ratios. 



this corresponds to a mass limit around 0.62M o for an age rang- 
ing between 50 and 90Myr (IChe et al.ll20TTl) . 



5.3. Limitations and possible improvements 

The error bars on the closure phase estimated by pndrs take 
into account two main contributions: (i) the statistical disper- 
sion of the data within the 100 scans in each individual file, 
and (ii) the fluctuations of the transfer function estimated from 
the calibration measurements. Up to now, the latter has gener- 
ally been dominant in the error bar budget, so that our obser- 
vations are not expected to be photon-noise limited. However, 
based on the PIONIER data obtained so far, we note that the 
sensitivity to off-axis companions generally improves when in- 
creasing the integration time and when observing brighter tar- 
gets, which are two characteristics of photon-noise limited ob- 



servations. Determining whether our observations are actually 
photon-noise limited or not is beyond the scope of this paper. 

One crucial aspect for high-precision interferometric obser- 
vations in general is the calibration. Although we have not noted 
suspicious variations in the closure phase tranfer function esti- 
mated from the calibration stars (as long as the weather condi- 
tions were decent), there might be some systematic effects when 
the calibration stars are located too far from the science target, 
especially at high air mass. Another possible question is the pres- 
ence of unknown faint companions around the chosen calibra- 
tors. Otherwise, calibration is usually not an issue: the closure 
phases of all our calibrators are generally consistent with zero 
within the error bars, and also with the science target's closure 
phases in case of a non-detection. 

Our experience with the PIONIER data shows that four to 
five OBs of good quality obtained at different hours angles are 
required to get robust results with our search method. While two 
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Fig. 6. Same as in Fig.[5]for the case of Regulus. In addition to the immediate blind test results (diamonds), for which the companion 
is deemed to be detected only when the contrast and position of the global best fit are compatible with the input parameters of the 
fake companion, we introduce another detection criterion independent of the position of the global best fit (asterisks), which is more 
suited for short integrations (see text). 



or three OBs are generally enough to reach a reasonably good 
dynamic range, the positions of the detected companions fre- 
quently remain ambiguous due to insufficient u, v plane cover- 
age. Furthermore, accumulating several OBs reduces the effect 
of possible systematics such as drifts in the closure phase at 
the edges of the wave band (which we could see in some data 
sets). The typical observing sequence to perform the astrome- 
try of faint companions around bright unresolved stars therefore 
consists in five observations of the target interspersed with obser- 
vations of reference stars, using at least two different references 
to reduce the systematic effects related to calibration. A large- 
scale survey for companions can, however, rely on two to three 
calibrated observations, which nowadays takes about one hour 
with PIONIER. 

The ultimate limitation to the dynamic range of PIONIER 
cannot be predicted yet, because on one hand, the level of the 
systematic errors cannot be precisely estimated from our limited 
data set, and because, on the other, the instrument and data re- 
duction software are expected to undergo significant upgrades in 
the future. In particular, we are currently investigating the pos- 
sibility of scanning the fringes more rapidly in the destructive 
read-out mode of the PICNIC camera, in an attempt to better 
stabilise the fringes during each individual scan. 

6. Conclusion 

In this paper, we have described our method to search for 
faint companions around bright stars with the PIONIER instru- 
ment at VLTI. We have demonstrated that dynamic ranges up 
to 1:500 can be reached when performing deep integrations (a 
few hours, including calibration measurements) on mostly unre- 
solved targets, while dynamic ranges around 1:200 can be rou- 
tinely achieved with three OBs of good quality, which nowa- 
days takes about one hour. We validated our search method by 
performing a double-blind test study, where fake companions of 
various flux ratios were introduced in our data sets. 

Our observations of Fomalhaut, tau Cet, and Regulus re- 
vealed the absence of companion within the 100 mas search re- 



gion, with median dynamic ranges at 3<x respectively of about 
1:430, 1:290, and 1:180. At 90% completeness, the dynamic 
ranges become 1:360, 1:220, and 1:130. In the first two cases, 
these non-detections confirm that the near-infrared excess emis- 
sions of about 1% previously detected around these two stars 
with high-precision visibility measurements at the VLTI and 
CHARA arrays are related to an extended source of emission 
rather than a point-like source. For del Aqr, we obtained a di- 
rect detection of a faint companion (flux ratio of 2.05 x 10~ 2 + 
0. 16 x 10~ 2 ) located at about 40 mas from the primary star, which 
was already known from astrometric and radial velocity mea- 
surements but was poorly constrained up to now. 

Based on these results and on the foreseen near-term de- 
velopments of PIONIER, we expect that dynamic ranges up to 
1 : 1000 should be reachable in the future. Whether PIONIER will 
eventually be able to meet its most ambitious goal, i.e., reaching 
the required dynamic range (~1 :5000) to directly detect hot giant 
exoplanets around unresolved target stars, is not clear yet. 

Acknowledgements. PIONIER is funded by Universite Joseph Fourier (UJF, 
Grenoble-1) with the programme TUNES-SMING, the Institut de Planetologie 
et d'Astrophysique de Grenoble (IPAG, ex-LAOG), and the Institut National 
des Science de l'Univers (INSU) with the programmes "Programme National 
de Physique Stellaire" and "Programme National de Planetologie". PIONIER 
was developed by the CRISTAL instrumental team of IPAG in collaboration 
with R. Millan-Gabet (NExScI) and W. Traub (JPL). The authors want to 
warmly thank the VLTI team. We acknowledge support from the French National 
Research Agency (ANR) through project grant ANR10-BLANC0504-01. This 
work made use of the Smithsonian/NASA Astrophysics Data System (ADS) and 
of the Centre de Donnees astronomiques de Strasbourg (CDS). 

References 

Absil, O., di Folco, E., Merand, A., et al. 2008, A&A, 487, 1041 
Absil, O., Le Bouquin, J.-B., Lebreton, J., et al. 2010, A&A, 520, L2 
Absil, O. & Mawet, D. 2010, A&A Rev., 18, 317 
Absil, O., Mennesson, B., Le Bouquin, J.-B., et al. 2009, ApJ, 704, 150 
Baraffe, I., Chabrier, G., Allard, E, & Hauschildt, P. H. 1998, A&A, 337, 403 
Berger, J. P., Haguenauer, P., Kern, P., et al. 2001, A&A, 376, L31 
Boccaletti, A., Augereau, J.-C, Baudoz, P., Pantin, E., & Lagrange, A.-M. 2009, 
A&A, 495, 523 

Borde, P., Coude du Foresto, V., Chagnon, G., & Perrin, G. 2002, A&A, 393, 
183 



10 



Absil et al.: Searching for faint companions with VLTI/PIONIER 



Bouchet, P., Schmider, F. X., & Manfroid, J. 1991, A&AS, 91, 409 
Chauvin, G., Lagrange, A.-M., Dumas, C, et al. 2004, A&A, 425, L29 
Che, X., Monnier, J. D., Zhao, M., et al. 2011, ApJ, 732, 68 
Di Folco, E., Absil, O., Augereau, J.-C, et al. 2007, A&A, 475, 243 
Di Folco, E., Thevenin, F, Kervella, P., et al. 2004, A&A, 426, 601 
Duvert, G., Chelli, A., Malbet, F, & Kern, P. 2010, A&A, 509, A66 
Goldin, A. & Makarov, V. V. 2007, ApJS, 173, 137 

Hinkley, S., Carpenter, J. M., Ireland, M. J., & Kraus, A. L. 201 1, ApJ, 730, L21 

Kalas, P., Graham, J. R., Chiang, E., et al. 2008, Science, 322, 1345 

Kervella, P., Thevenin, F, Di Folco, E., & Segransan, D. 2004, A&A, 426, 297 

Lacour, S., Tuthill, P., Amico, P., et al. 201 1, A&A, in press 

Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, in Proc. SPIE, Vol. 7734, 

Optical and Infrared Interferometry II, ed. W. C. Danchi, F. Delplancke, & 

J. K. Rajagopal 

Lagrange, A.-M., Desort, M., Galland, F, Udry, S., & Mayor, M. 2009a, A&A, 
495, 335 

Lagrange, A.-M., Gratadour, D., Chauvin, G., et al. 2009b, A&A, 493, L21 
Le Bouquin, J.-B., Berger, J. P., Lazareff, B., et al. 201 1, ArXiv e-prints 
Marois, C, Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 
Mawet, D., Mennesson, B., Serabyn, E., Stapelfeldt, K., & Absil, O. 2011, ApJ, 
submitted 

McAlister, H. A., ten Brummelaar, T. A., Gies, D. R., et al. 2005, ApJ, 628, 439 
Mennesson, B., Serabyn, E., Hanot, C, et al. 201 1, ApJ, in press 
Merand, A., Borde, P., & Coude du Foresto, V. 2005, A&A, 433, 1155 
Monnier, J. D. 2003, in EAS Publications Series, Vol. 6, Observing with the 

VLTI, ed. G. Perrin & F. Malbet, 213-226 
Oppenheimer, B. R. & Hinkley, S. 2009, ARA&A, 47, 253 
Zhao, M., Monnier, J. D., Che, X., et al. 2011, PASP, 123, 964 



